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The plateau- non- plateau transition in the antiferromagnetic- (ferromagnetic) „ polymerized S = 
1/2 XXZ chains under the magnetic field is investigated. The universality class of this transition 
belongs to the Brezinskii-Kosterlitz-Thouless (BKT) type. The critical points are determined by 
level spectroscopy analysis of the numerical diagonalization data for 4 < p < 13 where p(= n + l) 
is the size of a unit cell. It is found that the critical strength of ferromagnetic coupling decreases 
with p for small p but increases for larger enough p. It is also found that the plateau for large 
p is wide enough for moderate values of exchange coupling so that it should be easily observed 
experimentally. This is in contrast to the plateaus for p — 3 chains which are narrow for a wide 
range of exchange coupling even away from the critical point. 

KEYWORDS: magnetization plateau, level spectroscopy, exact diagonalization, BKT transition, twisted boundary 
condition 



§1. Introduction 

Magnetization plateaus in one dimensional Heisenberg 
chains are attracting much attention as the field induced 
spin gap states. On the plateau, some spins are partly 
quenched by the magnolia fir^ld_wJdk iJiajcmaining spins 
form a spin gap stateJD&WliBBB&fciBBfca> It is evi- 
dent that such magnetization plateaus can be realized in 
various non-homogeneous Heisenberg chains. Oshikawa, 
Yamanaka and Afflecku^ proposed the necessary condi- 
tion for the magnetization plateaus as p(S — m z ) = q = 
integer where p is the periodicity of the magnetic ground 
state in the thermodynamic limit, S is the magnitude of 
the spin, and m z is the magnetization per site. 

In this context, the spin chains consisting of pe- 
riodic arrays of ferromagnetic and antiferromagnetic 
bonds have been widely investigated. The S = 1/2 
ferromagnetic-antiferromagjietic alternating chain has a 
spin gap ground stateOJ'LLa'EJf for arbitrary values of 
coupling constants. This state can be regarded as the 
plateau state with zero magnetization (p = 2, S — 
1/2, q = 1 and m z = 0). In the isotropic case, this 
energy gap survives as the well-known Haldane gap even 
if the strength of the ferromagnetic bonds j is infinite.t!-f 

There appears a magnetization plateau in the magneti- 
zation curve of the S = 1/2 ferromagnetic- ferromagnetic- 
antiferromagnetic trimerized chain at 1/3 of the satura- 
tion magnetization (p = 3, S = 1/2, q = r l and m z = 
1/6) .cP In this case, Kitazawa and Okamatca' have shown 
that the plateau vanishes by the Brezinskii-Kosterlitz- 
Thouless (BKT) type transition for j > j c ~ 15.4 where 
the energy unit is the strength of the antiferromagnetic 
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bonds. 

In the present work, we investigate the plateau-non- 
plateau transition for general values of p using the 
method of ref. [)| Intuitively, the quantum fluctuations 
are expected to be suppressed as the number of the fer- 
romagnetic bonds increases. Therefore it is natural to 
expect that the magnetization plateau, which is an es- 
sentially quantum phenomenon, becomes less favourable 
as the number of the ferromagnetic bonds increases. Sur- 
prisingly, our numerical calculation shows that this is in- 
correct. We discuss the physical origin of this unexpected 
behavior. 

This paper is organized as follows. In the next section, 
the model Hamiltonian is defined and the numerical re- 
sults of the magnetization curve and the phase diagram 
of the plateau-non-plateau transition is presented. The 
physical reason of the enhancement of the plateau state 
for large p is also explained in this section. The final 
section is devoted to a summary and discussion. 

§2. Numerical Results 

2.1 Model Hamiltonian 

The Hamiltonian of the antiferromagnetic- (ferromagnetic)„ 
polymerized S = 1/2 XXZ chain in the magnetic field is 
given by 

L n-l L 

K = ~ 3 ^ ^ + ^ ft{n+l)l+n(A) 

1=1 i=0 1=1 

(n+X)L 

-gp B H S i- C 2 - 1 ) 
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where 



Hi(A)=SfS, 
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The ferromagnetic coupling constant, number of the fer- 
romagnetic bond and anisotropy parameter of the anti- 
ferromagnetic bond are represented by j, n and A, re- 
spectively. The ferromagnetic bonds are assumed to be 
isotropic. The magnetic field, the electronic g-factor and 
Bohr magneton are denoted by H, g and /ib > respectively. 
This model has spatial periodicity p = n + 1. In the fol- 
lowing, we take gfj,B = 1 as unit. 

In the limiting case of j — -> oo, this model tends to 
the S = p/2 uniform antiferromagnetic XXZ chain, so 
that no plateau is expected. Therefore the BKT type 
plateau-non-plateau transition is expected to take place 
at a finite critical value of j. 



numerical diagonalization data for general values of p. 
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2.2 Magnetization plateaus 

The schematic magnetization curve of the present 
model is shown in Fig. 1. In the present paper, we 
concentrate on the highest plateau with q = 1 and 
m z = 1/2 — 1/p. The total magnetization M = m z 
on the plateau is given by M = M p = pi (1/2 — 1/p) = 
M s (l — 2/p) where M s = pL/2 is the saturation mag- 
netization. In Fig. 1 the critical fields H c q corresponds 



M/M 



1-2/p 




Fig. 1. Schematic magnetization curve of the present model. 



to the energy gap in the absence of the magnetic field, 
H c i, to the energy gap between the lowest energy with 
M z = M p and that with M z = M p — 1, H e2 , the energy 
gap between the lowest energy with M z = M p + 1 and 



that with M z = M p . Finally H s is the saturation field. 
The lower part of the magnetization curve is represented 
by the dashed line where lower plateaus may appear. 

We use an exact diagonalization method to calcu- 
late the magnetization curve. As an example, the j- 
dependence of the 3 critical fields H c q, H c i and H C 2 for 
p = 4 and A = 1 are shown in Fig. 2 by the symbols 
and o, respectively, which are obtained by the Shanks 'tZP 
extrapolation of the results of the Lanczos exact diago- 
nalization for AL = 8, 16 and 24 with periodic boundary 
condition. The plateau seems to close at a finite value 
of j. Because the non-plateau state corresponds to the 
gapless spin liquid and the plateau state to the spin gap 
state, this transition is expected to be the BKT-type 
transition. In the next subsection, we determine the 
critical point of the plateau-non-plateau transition from 



Fig. 2. The critical fields for p = 4 and A = 1. The critical fields 
H c q,H c i and H c 2 are represented by •, O and o. 



2. 3 BKT critical point 

It is difficult to estimate the BKT critical point 
precisely from standard finite size analysis of the nu- 
merical calculation data. To circumvent this diffi- 
culty Nomura and KitazawaH proposed to use level 
spectroscopy methodE3' with twisted boundary con- 
dition. This method has been successfully applied 
to the ferromagnetic-ferromagnetic-antiferromagnetiij 
trimerized Heisenberg chain by Kitazawa and Okamoto.Er 
In the present work, we employ this method for the 
present model. Here we do not explain how and why 
this method works, because methodological details are 
well described in ref. [| 

The finite size critical point is determined from the 
crossing point of AE 0: 2 with the lower of Ai?^ 2 or 
A£^/ 2 defined by, 

AE , 2 = i{£ (M p + 2,0,1) + E (M p - 2,0,1)} 

-£ (M P ,0,1). (2.3) 



and 



AE l/2,0 - E 



TBC 



(Mp,l)-£(M P ,0,1) 



/2,0 



= E 



TBC 



(Af p ,-1)-E(M P ,0,1) 



(2.4) 



(2.5) 



where Eo(M z , k,P) is the lowest energy under a periodic 
boundary condition with magnetization M z , wave num- 
ber k and parity P. The energy E TBC (M z ,P) is the 
lowest energy with the twist boundary conditions with 
magnetization M z and parity P. 

To confirm the reliability of this method the following 
average 

x i/2.o(L) + 3x1, 2 (L) 



(2.6) 

is close to 0.5 at the critical pointsi-' where x c 1 J 2a (L) 
are the scaling exponents corresponding to E C ^J 2 and 
defined by 



'1/2,0 



2^ A£ V2,o- 



(2.7) 
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Here v s is the spin wave velocity given by 

v s = lim ±-[E MpM (L)-E M (L)], 

L—>oc Z7T 



(2.8) 



where EM p ,ki_{L) is the energy of the excited state with 
wave number k\ = 2jS and M z = M p . The results are 
shown in Fig. 3 for p = 4 which confirms that x = 0.5 
holds with good accuracy. 

Figure 4 shows the behavior of A£^o,2 and AEf, 2 for 
L = 6 unit cells and p = 4. From the crossing point, 
we obtain j c (L) = 6.310 for L = 6 unit cells. The BKT 
transition point for the infinte system can be obtained 
by extrapolating from L — 2, 4 and 6 to N — > oo as 
jc = 6.260 assuming the extrapolation formula j c (L) — 
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Fig. 3. The extrapolated values of the averaged scaling dimension 
x on the critical points for p = 4. The solid line is x = 0.5 
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Fig. 4. The j'-dependence of the energies A_Eo,2 and AE^ 2 Q rep- 
resented by o and •, respectively, for p = 4, A = 1 and L = 6. 



The phase boundaries are shown in Fig. 5 for p = 3, 4 
and 5 by □, • and o, respectively. Extrapolation is made 
from L = 2, 4 and 6 for p = 4 and 5. The results for p = 3 
are taken from r ef. § Fr om this figure, we see that the 
critical points j c does not decrease monotonically with 
an increase of periodicity p. Therefore we further calcu- 
late the cases of longer periodicity up to p = 13. The p 
dependence of j c is shown in Fig. 6 for A = 1.0,0.5,0.0 
and —0.5. The different symbols represent the critical 
points for different system size L. The data for p = 2 
is also included for A = ,=.0.5. When A > and p = 2 
then j c equals infinitylL§ \3 For p — 4 and 5, the calcu- 
lations are carried out for L = 2,4 and 6 and for p = 6 
and 7, L = 2 and 4. The extrapolation to L — > oo is 
made for these cases. Nevertheless, the L-dependence of 
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Fig. 5. The phase boundaries on the A — j'-plane for p = 3, 4 and 
5 represented by □ • and o , respectively. The lines are fitted. 



Jc 

15 
10 






1 












A 








. L=2 






A 








o L=4 






+ 








a L=6 














- L=8 














No Plateau 














A=1.0 

• 


• 


• 

• 




• 






• 




• ■ 










+ * A=0.5 


• 


• 










• 
















=0.0 


• 








■f 


• 

• 


• 




+ 


+ 




+ 




A=- 


0.5 
• • 


+ 

-f 


+ 




♦ 


, . . • ' 

Plateau 





12 P 



Fig. 6. The p-dependence of the critical point j c for A = 
1.0,0.5,0.0 and —0.5 from top to bottom. The different sym- 
bols distinguish j c 's for different system sizes L. The + symbols 
are the values of ^^extrapolated to L — ► oo. The repq^jfe^or p = 3 
is taken from ref. E* For p = 2, j c — > co for A > 0.1 



the critical point becomes weaker as the periodicity p be- 
comes larger as shown in Fig. 6. Therefore for p > 8, we 
can safely estimate approximate values of the thermody- 
namic critical point j c from that for L = 2. Surprisingly, 
the p-dependence of the critical point is strongly non- 
monotonic. The decrease of j c with p for small values 
of p is consistent with intuition that the increasing the 
number of the ferromagnetic bonds suppresses the quan- 
tum fluctuation. The magnetization plateau, which is 
the quantum origin, becomes less favoured. 

With the further increase of p, however, j c increases 
again almost linearly with p. This phenomenon can be 
explained by the following physical picture. On the mag- 
netization plateau M z = M p , the singlet pairs are lo- 
cated on the antiferromagnetic bonds and other spins 
are alligned in the direction of the magnetic field. This 
is schematically shown in Fig. 7(a). If the magne- 
tization becomes M z = M p + 1, one singlet pair is 
broken as shown in Fig. 7(b). The energy difference 
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Eo(M p + 1) — Eo(Mp) determines the upper edge of 
the plateau H C 2- Because the additional up spin is lo- 
calized around the antiferromagnetic bonds, the criti- 
cal field H C 2 tends to a constant value of the order of 
the antiferromagnetic coupling constant ~ 0(1) as p in- 
creases. On the other hand, if the magnetization be- 
comes M z = M p — 1 , one of p — 2 spins in the ferromag- 
netic segments is inverted as shown in Fig. 7(c). In the 
absence of coupling to the 2 boundary spins forming a 
local singlet, the segment of the remaining p — 2 spins, 
which contain one inverted spin, form a state with total 
spin S = and magnetization S z = 1 -^ L — 1. Since 
the ferromagnetic bonds are isotropic this state has the 
same energy as the fully polarized state with S = 
and S z = 2^2. The energy change E (M p ) - E (M p - 1) 
is due to the interaction of the inverted spin with the 
singlet pair on the antiferroamgnetic bonds at the both 
ends of the ferromagnetic segment. The probability 
that the inverted spin stays on the boundary sites of 
the segment is proportional to — i^. Hence, the energy 
diference Eq(M p ) — Eq(M p — 1), which determines the 
lower edge of the plateau H c \ is estimated as O(j^) for 
large enough p. Because the plateau-non-plateau tran- 
sition takes place around the point H c \ ~ H C 2, we find 
j c oc p — 2 for large enough p. This explains why j c in- 
creases linearly with p for large p. To substantiate this 
argument, H c \ and H C 2 are calculated numerically by 
exact diagonalization for L = 2 and A = 1.0. Figure 8 
shows that H c \ scales with for p = 10, 11, 12 and 13 
from top to bottom although there is some deviation as 
the critical point is approached. Figure 9 shows H C 2 is 
approximately independent of p and tends to a constant 
value for large j. These results show that the increase of 
the critical point for large p is a consequence of delocal- 
ization of the down spin in a wide region of length p. In 
Fig. 10, 18 plotted against p. It is evident that 
tends to a constant value for large enough p. 
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Fig. 7. The spin configurations (a) on the magnetization plateau 
M z = M p , (b) at the magnetization M z 



1 and (c) at 



the magnetization M z = 
ovals form local singlets. 
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§3. Summary and Discussion 

The magnetization plateaus at m z = 1/2 — 1/p in 
the S — 1/2 antiferromagnetic-(ferromagnetic)„ poly- 
merized XXZ chains are investigated by exact diagonal- 
ization of finite size systems. The plateau-non-plateau 
BKT transition points are determined precisely by the 
leveLsxysctroscopy method with twisted boundary condi- 
tionM 

The p-dependence of the critical point is not a mono- 
tonic function. For small p, the critical point decreases 
with p in agreement with the simple-minded intuition 
that an increase of the number of the ferromagnetic 
bonds suppresses the quantum fluctuation and the mag- 
netization plateau, which is the quantum origin, becomes 
less favoured. With a further increase of p, however, j c 
increases again. This can be explained by the delocaliza- 
tion of down spin within the ferromagnetic segment with 
length p — 2. 

The plateau for p — 3 is narrow for a wide range of 



Fig. 8. The ^ dependence of H cl for p = 10, 11, 12 and 13 
from top to bottom, H c \ scales as . 



3 < jcO 1 This is the reason why so far this simplest 
example of the magnetization plateau has not even ob- 
served experimentally. Actually, no visible plateau is ob- 
served in the magnetization curve of 3CuCl2- 2dx which 
has j ~ 4.5 << j c ErEf In contrast, the magnetization 
plateaus for large p is fairly wide compared to that for 
p = 3 for moderate values of j. As an example, we show 
the magnetization curve for p = 12, A = 1 and j = 4.0 
obtained by exact diagonalization for L = 2 in Fig. 11. 
It is evident that this plateau is wide enough for experi- 
mental measurement, although the value of j c for p = 12 
is less than that for p = 3. Thus a systematic synthesis 
of materials with broad spatial periodicity is hoped to 
experimentally confirm our predictions. 

Chen and Hida thank K. Okamoto for fruitful discus- 
sion. We also thank A. Kitazawa and K. Okamoto for 
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Fig. 9. The j'-dependence of H c2 for p = 10, 11, 12 and 13. H c2 
is almost independent of p and the tends to a constant value for 
large j. 
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supplying us the numerical data of ref. ^. The numeri- 
cal calculation is performed using the program package 
TITPACK version 2 coded by H. Nishimori on HITAC 
S820 and SR2201 at the Information Processing Cen- 
ter of Saitama University and the FACOM VPP500 at 
the Supercomputer Center of Institute for Solid State 
Physics, University of Tokyo. This work is supported by 
a Grant-in- Aid for Scientific Research from the Ministry 
of Education, Science, Sports and Culture of Japan and a 
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Fig. 10. The p-dependence of j c /(p — 2). 




Fig. 11. The magnetization curve for p = 12, A = 1 and j = 4 
obtained by exact diagonalization with L = 2. 



